Upper limits on electric dipole moments 
of r-lepton, heavy quarks, and VF-boson 
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Abstract 

We discuss upper limits on the electric dipole moments (EDM) of the r-lepton, heavy 
quarks, and VF-boson, which follow from the precision measurements of the electron and 
neutron EDM. 



1 Introduction 

Strict upper limits on the electric dipole moments of common elementary particles, electron 
and proton, were derived from spectroscopic, almost table-top experiments [HISIE]. As to the 
neutron EDM, the best upper limit on it was obtained as a result of reactor experiments lasting 
many years [4| (they say that the searches for the neutron EDM killed more theories than any 
other experiment in the history of physics). And at last, the result for the muon EDM follows 
from the measurements at the dedicated muon storage ring [5]. These results are summarized 
in Table 1. 
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(0.7 ±0.7) X 10-27 ^ij 


< 0.8 X 10-24 ^2, 3 


< 0.29 X 10-25 4] 


(0.37 ±0.34) X 10-18 [5 



Table 1 

As to the dipole moments of the r-lepton and heavy quarks, upper limits on them have 
been obtained up to now from the analysis of high-energy experiments. 

The approach pursued here is based on the precision results [HH]. We establish upper limits 
on the dipole moments of the r-lepton and heavy quarks through the analysis of their possible 
contributions to the electron EDM. This is a clean theoretical problem for the r-lepton. For 
heavy quarks and VT-boson these our results are of rather qualitative nature. Additional upper 
limits on the dipole moments of heavy quarks, and W^-boson, also qualitative ones, are derived 
by the analysis of their possible contributions to the neutron EDM. 
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2 Dipole moment of r-lepton and electron EDM 
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Figure 1 



We start with the analysis of the contribution of the r EDM d^- to the electron dipole 
moment de- This contribution is described by the diagrams of the type presented in Figs. 
la,b,c,d. Here the loop is formed by the r line, and the lower solid line is the electron one. 
The upper wavy line corresponds to the external electric field. The crossed vertices refer to the 
electromagnetic interaction of the r EDM 
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Of course, all six permutations of the electromagnetic vertices on the electron line should be 
considered. The contributions of diagrams lb and Ic are equal. 

This problem is similar to that of the contribution by the light-by-light scattering via muon 
loop to the electron magnetic moment [6]. The general structure of the resulting contribution 
to the electron EDM is rather obvious (to the leading order in me/m,-): 
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where a is some numerical factor (hopefully, on the order of unity). The factor me originates 
from the necessary helicity-flip on the electron line; then l/irij. is dictated by dimensional 
arguments. 

Diagram la corresponds to the matrix element <e\fa^^T\e> = Cua^'^u. We use dimensional 
regularization with d = A — 2e dimensions and the method of regions (see the textbook [7]). 
Only the region where all three loops are hard (loop momenta ~ m^) contributes to the leading 
power term ([2]); therefore, there are no logarithms ln(mT-/me) (contributions of regions with 1 
or 2 hard loops are suppressed by an extra factor (mg/m^)^). In this hard region, the problem 
reduces to 3-loop vacuum integrals with a single mass belonging to the simpler topology 
Bm [8]. We perform the calculation in arbitrary covariant gauge, and use the Reduce package 
Recursor [8j to reduce scalar integrals to two master integrals. Gauge-dependent terms 
cancel, and we get 
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and ( is the Riemann ^-function. All divergences cancel, and we arrive at the finite contribution 
to a: 

a. = |c(3)-^. (5) 

In order to calculate the contribution of Fig. lb,c,d, we expand the corresponding initial 
expressions in the external photon momentum q up to the linear term. The EDM vertex 
contains £i^'"^f^-^ we put this factor aside, and calculate tensor diagrams with four indices. After 
summing all diagrams, the result is finite; now we can set e ^ 0, and multiply by e'""''^ (cf. [H]). 
The result has the structure of a tree diagram with the electron EDM vertex s^"'^^ qya^p- The 
gauge- dependent terms in it cancel (exactly in d), as well as the divergences. This contribution 
to a is 

a2 = ^C(3)-l. (6) 

As an additional check of our programs, we have reproduced the leading power term in the 
contribution to the electron magnetic moment originating from the light-by-light scattering via 
the muon loop (formula (4) in [6]). 

The final result for the numerical coefficient is 

15 31 

a = ai + as = - C(3) - - = 1.924 . (7) 

With this value of a, the discussed contribution to the electron EDM is 

= 6.9 X lO^^^d^. (8) 

Combining this result with the experimental one [1] (see Table 1) for the electron EDM, we 
arrive at 

d^/e = (1± 1) X lO^^^cm. (9) 

In fact, the results ([5]) and (I6l) refer to somewhat different regions of incoming momenta. 
For ([6]) all the three momenta are hard, on the order of magnitude about m^-, but for ([5]) only 
two of them belong to this region, and the third one, that of the outer photon, is soft, of 
vanishing momentum. Still, one may expect that the effective EDM interaction is formed at 
momenta much higher than m,-, so that this difference is not of much importance. Besides, the 
contribution of diagram la is anyway numerically small. Thus, result ([5]) is valid at least for 
all momenta about m-r ~ 1 - 2 GeV. 

The upper limits on the r EDM derived from the accelerator experiments [TOt [TTl [T2l [13] 
belong to the interval of 10^^^ — 10^^'' e-cm, so that our result ([9]) formally does not improve 
them. However, all those accelerator data refer to much larger typical momenta of the photon, 
from 10 to 200 GeV. 



3 Dipole moments of heavy quarks and electron EDM 

The information on the dipole moments dq of heavy quarks can be obtained from the diagrams 
analogous to those presented in Figures 1, but with quarks in the fermion loops, instead of 
the r-lepton. We do not consider the gluon corrections to the quark loops, and confine to the 
simple estimate for them, following from the analogy between the two problems (see ([2])): 




here rriq is the quark mass, Q is the quark charge in the units of e {Qb = —1/3, Qc,t = 2/3). 
As to the overall numerical factor here, we put ~ 1. The corresponding estimates are 
straightforward. Here and below we assume the following values for the quark masses: 

TUc = 1.25 GeV; = 4.5 GeV; mt = 175 GeV. 

Thus obtained upper limits are: 

dje < 3 X 10-^'' cm; 4/e < 7 x 10"^^ cm; dt/e < 4 x 10"^^ cm . (11) 

4 Dipole moments of heavy quarks and neutron EDM 

The dipole moment of a heavy quark Q generates the EDM of a light quark q via the following 
diagram [I4j: 




The scattering amplitudes and corresponding contributions to the dipole moments of light 
quarks are: 

Ml = -r^TTT I^Q^I' 6^ ^ - 2^ ^Q^75a,.« F,^ , (12) 
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for 



Q = t, q = d. 



Here Vq^ are the corresponding coefficients of the Kobayashi-Maskawa matrix; 6'^ is the Wein- 
berg angle, sin^^^ = 0.23. The structures of relations f|T3l) and f|T5l) are different since 
^c,fe ^ ^VK; while > In fact, to simplify final expressions, we confine to the leading 

order in 1/m^ in ([T^ . ([T^ . and even assume that nit ^ in (fT^ - ([IS])- In both cases 



the diagrams are logarithmically divergent, so that we have to introduce a cut-off A at high 
momenta. 

This divergence is caused by the term k^kiy/m'^ in the numerator of the W Green function 
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Therefore, relations (fT51) and (fT^ differ 



The mentioned term k^ky/m^ was omitted in 
essentially from the corresponding results of [Hj. 

In our estimates here, we assume that both /m^) — 2 and /m^) — 2 are on the 
order of unity. As to the light quark dipole moments, we assume (in the spirit of the constituent 
quark model) that they are on the same order of magnitude as the neutron EDM: 
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Of course, both these assumptions make the corresponding estimates less definite than those 
based on electron EDM. The results of our estimates are presented in the last line of Table 2. 
In its previous line we repeat for comparison the data already given in ( fTTi) . In both sets of 
estimates the strong interaction of quarks is neglected. 
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dq/e, cm 
derived with de 


< 3 X 10-16 


< 7 X 10-15 


< 4 X 10-14 


dq/e, cm 
derived with dn 


^ 10-15 


^ 2 X 10-12 


^ 5 X 10-15 



Table 2. Limits on quark dipole moments from electron and neutron EDMs 

It is only natural that for the most heavy t-quark the bound from dn dominates. 

The bounds on the dipole moments of c- and 6-quarks derived from the high-energy experi- 
ments p!3| [T5] are on the order of 10~i^ e-cm, so that they are more strict than the limits given 
in Table 2. However, as noted above already, the high-energy data refer to different, much 
larger typical momenta of the photon. 

In conclusion of this section, we mention the investigation of the radiative decay b sj a.s a. 
probe of the t -quark EDM [16] . The discussed contribution to the decay amplitude is described 
by the diagram in Fig. 3. This diagram is also logarithmically divergent. Its contribution can 
be easily obtained as follows. In the limit — > rris, the structure of discussed matrix element 
should coincide with that of ([Hj) (up to the obvious changes — >• m^g) and \Vtd\'^ V-i*gVtb)- 
On the other hand, in the limit —>■ 0, the outgoing s-quark should be left-handed. Therefore, 
the matrix element of the b ^ sj decay (with the internal t-quark) should look as follows: 
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Figure 3 



a ^^^^^ mbrrit ( ^ 



64vr sin^ ^ ( - 2 ) ^.(1 " 75) F,. . (16) 



The result presented in [l6j is 

df/e < 10"^^ cm. (17) 

5 Electric dipole moment of Vl^-boson 




eve 
Figure 4 



One more contribution to the electron and neutron dipole moments can be given by the 
EDM dw of H^-boson. This effect was pointed out and investigated long ago [TTl [18] . 

It is convenient to start the discussion with the electron EDM. The effect is described by 
diagram presented in Fig. 4. Here as well the crossed vertex refers to the electromagnetic 
interaction of the W^-boson EDM (it is obvious from the diagram that here it is the dipole 
moment of W~\ In this case the EDM interaction is described by the Lagrangian 



' edm 



2mwidwFapWlWp . 



The corresponding matrix element is 
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With straightforward, though rather tedious calculations (somewhat simplified by employing 
the density matrix of polarized fermion), one arrives at the following result for the contribution 
of ly-boson EDM to electron dipole moment: 
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Here A is the cut-off parameter for the logarithmically divergent integral over virtual momenta 
in the loop. Putting (perhaps, quite conservatively) ln(A^/m^) ~ 1, one obtains with the 
experimental upper limit on the electron EDM [Ij (see Table 1), the following bound on the 
dipole moment of ly-boson: 



dw/e ^ 2 X 10 cm. 



(21) 



In the case of the W^-boson contribution to the neutron EDM, our line of reasoning somewhat 
differs from that of [Ul US]. We note first of all that the electron mass me does not enter 
explicitly matrix element (fT9|) . It arises in the result (120|) only as the mass of an external 
fermion, via the Dirac equation pu = rrieU. Therefore, there are all the reasons to expect that 
the contribution of dw to the neutron EDM will be proportional to the neutron mass m„, i.e. 
enhanced as compared to (|20l) by three orders of magnitude. In this case, the forward scattering 
amplitude of the virtual W^-boson can be written in a general form as follow^: 



Mn(p)7/.(1 +75) kg{k'^) + ph{k^) 7^(1 + 75)m„(p) ; 



(22) 



here k is the total momentum of intermediate hadronic states. Of course, the invariant functions 
g and h depend in fact not only on /c^, but on (kp) as well. However, in our case fc^ ~ ^ 
(kp) ~ ninmiY, so that the dependence on (kp) can be safely ignored. By the analogous 
reason, in the usual case of the deep inelastic neutrino scattering, the structure with p in the 
corresponding amplitude is also omitted. In the present case, however, we should keep in 
amplitude (l22l) p, in addition to the common k, since after integrating over d'^k both structures 
give comparable contributions to the result. 

At last, the usual dimensional and scaling arguments dictate that asymptotically, for k"^ ~ 
m^, both functions g and h behave as follows: 
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.(P)^g, Hk^)^'^,. 

In particular, one can neglect the gluon corrections in these functions. Without any additional 
parameters, it is natural to assume that go, ho ^ 1. 

Now, the same calculations as those in the case of electron EDM, result in the following 
expression for the discussed contribution to the neutron dipole moment: 
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'^Compare with the corresponding structure Me(p)7/i(l+75) {k/k"^) 7iy (l+75)ue(p) in formula (|19p for electron. 



Then, with the result of [1] for the neutron EDM (see Table 1), we arrive at the following quite 
strict upper limit on the l^-boson dipole moment: 

c^Hz/e < 2 X 10^2^ cm . (24) 
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